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Abstract: Several studies have characterized the Yungas as a separate biogeographic unit, mainly based in floristic 
components. However, these characterizations were mainly qualitative and did not include faunal groups. The Yungas 
have been assumed as a region with rich floral and faunal diversity, but without testing how well they are described by 
animal distributions. Our study consists of a formal analysis of endemism based on distribution of small mammals in the 
southernmost portion of the Yungas. This area is biogeographically very interesting because the Yungas are comprised of 
discontinuous fragments of forests that extend into temperate arid and semiarid habitats. As a first approximation, we 
contrasted a group of volant species (bats) versus a group of non-volant species (marsupials). Our results show that small 
mammals are efficient indicators of endemism in Yungas. Eighty percent of the species of small marsupials included in 
the analysis supported the identified areas as being zones of endemism. Regarding bats, almost 55 percent of the species 
supported a designation of endemism. The results also show that the areas we considered are congruent with the botanical 
definition of the Yungas of northwestern Argentina and their subdivisions, an assumption that had not been previously 
tested with a formal quantitative method. We also found that non-volant species are better indicators of endemism than 
volant ones at regional scales, but volant species are better indicators than was previously thought. 
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It is well known that mammals are important as 
biogeographic indicators and their species, or higher 
taxonomic levels (genera, families, and orders), have 
been important in defining such areas at different scales. 
The Biogeographic Regions of the world have endemic 
orders, families, and subfamilies, while Subregions and 
minor divisions have endemic genera and species. 
Koopman (1976, 1981) and Proches (2005) found that 
bats are important for defining faunal regions, while 
other authors (Willig & Mares, 1989; Corbet & Hill, 
1992) suggested that bats are not useful for defining 
biogeographic zones or their limits within lowland 
tropical forest. Proches (2005), in his analysis of bat 
distributions at global scales, obtained regional 
groupings largely consistent with the zoogeographic 
regionalization model and the phytogeographic data are 
consistent with the zoological data. His results are 
coincident with that of Koopman (1976, 1981), who had 
already considered endemics among bats as defining 
faunal provinces of South America. On the other hand, 
Corbet & Hill (1992) pointed out that the distributions of 
non-flying mammals are more tied to the regional limits 
than are the distributions of flying mammals. They 
observed than the Australian and Asian faunas of volant 
species are more mixed than the terrestrial mammals and 
any limit is arbitrary. Similarly, Willig & Mares (1989) 
obtained little congruence between phytogeographic 
zones and distributional limits of bat species in 
Venezuela, so they argued that the distribution of bats 
were not limited to particular phytogeographic zones; 
moreover, they argued that no associations of bats can be 
used to define biogeographic zones. 

Descriptions of the Yungas have been essentially 
qualitative (Cabrera & Willink, 1973; Cabrera, 1976; 
Vervoorst, 1979) and based on floristic components (see 
Ribichich, 2002). For animals the Yungas were then 
assumed to be a natural region, but they did not test the 
congruence of species distributions of animals or plants 
with the putative regions. Moreover, the large latitudinal 
extension and the great altitudinal development of the 
mountains of northwestern Argentina region provide 
great environmental heterogeneity, leading several 
authors (e.g. Cabrera & Willink, 1973; Brown, 1995a, b; 
Hueck, 1978; Ojeda, 1999) to consider it as a non- 
homogenous unit, and to assume that its fauna was 
composed by species from the surrounding regions. 
However, other authors suggest that in the Yungas there 
are taxa that are more typical and common there than in 
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marginal areas, although they may not be exclusive 
(Cabrera, 1976; Morrone, 2000, 2001). 

The distribution of animals and plants in a region is 
the result of a complex interplay of historical, 
evolutionary, and ecological factors acting over different 
time and spatial scales (Hortal et al, 2010; Levin, 1992). 
A number of species with different dispersal capabilities 
which have similar distributional areas share a common 
history, besides similar ecological requirements. To 
share a common history implies that a number of taxa 
have been affected by the same historical processes that 
have operated for millions of years and that they have 
responded in similar ways. The congruent distributional 
areas of these taxa form an area of endemism. According 
to Platnick (1991) an area of endemism is defined by the 
areas of consistent distribution of two or more species 
spatially restricted to it. This definition implies sympatry, 
although not necessarily an exact congruence of the 
distributional areas of the species. Various delineations 
of areas of endemism over the years (Harold & Moi, 
1994; Morrone, 1994; Linder, 2001; Szumik et al, 2002) 
are based on the definition established by Platnick (1991) 
(Szumik $ Roig-Juñent, 2005), but only Szumik et al 
(2002) formally compared the degree of superposition 
between the distributional areas of species with the 
application of an optimality criterion. This comparison is 
fundamental considering that sympatry is the basic 
requirement to establish areas of endemism (Szumik & 
Roig-Juñent, 2005). 

The biogeographic assessment of a region is a 
starting point to formulate strategies for conservation of 
biodiversity. Consequently, in order to propose and 
develop programs of conservation and management, it is 
important to have quantitative descriptions of the 
different components and levels of biodiversity of a 
particular region (Ojeda, 1999). In this sense the areas of 
endemism are crucial (Vane-Wright et al, 1991; Linder, 
1995; Peterson and Watson, 1998). Previous studies 
dealing with identification of areas of importance for the 
conservation of mammals in Northwestern Argentina 
(NWA) that use quantitative methods over data of 
occurrence do exist (Ojeda et al, 2003; Tabeni et al, 
2004). Also, the areas of endemism are considered the 
fundamental units in the study of historical biogeography, 
because they are the entities to be compared in the 
biogeographic analyses. Finally, the identification of 
areas of endemism can be useful as additional evidence 
for studies that involve species of uncertain taxonomic 
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status, and to predict the occurrence of species in areas 
where they have not been recorded (Díaz Gómez, 2007). 

We have applied an explicit method to mammal 
distributions based on grids in order to identify areas of 
endemism. As a first approximation, we have contrasted 
selected small mammals, volant (bats) versus a group of 
non-volant (marsupials), from the Andean Subtropical 
Forests of NWA. We have evaluated whether small 
mammals in general are useful indicators of areas of 
endemism, whether the areas recovered are congruent 
with perceived natural areas (as botanically defined), and 
whether non-volant species are better indicators than 
volant one at the regional (meso) scale. 


1 Materials and Methods 


1.1 Study area 

The study area (Figs. 1, 2) includes the provinces of 
Jujuy, Salta, Tucuman, Catamarca, and Santiago del 
Estero in northern Argentina. This region is a complex 


Catamarca 





mosaic of habitats that intermingle, where large 
phytogeographic units converge, and some reach their 
boundaries of occurrence (Cabrera, 1976). The Yungas 
phytogeographic province (Cabrera & Willink, 1973; 
Cabrera, 1976) forms a narrow slash from Venezuela 
south to NWA, along more than 4,000 km, from 9° north 
latitude, and gradually disappearing at 28” south latitude 
in the province of Catamarca (Cabrera 8 Willink, 1973; 
Hueck, 1978). The southernmost part of the Yungas 
forests is recognized as a distinctive biogeographic and 
ecological unit (Cabrera, 1976). Moreover, Ojeda et al 
(2008) emphasized the distinctive character of this 
fraction of the Yungas, located in Argentina between 22° 
and 28° south latitude (from the border line with Bolivia, 
at 22° south latitude, to the centre of Catamarca; Cabrera, 
1976), defining it as a “peninsula” because of its position 
in relation to the habitat matrix that surrounds it. This 
portion of the Yungas is then unique and biogeogra- 
phically very interesting. The entire zone consists of an 


Fig. 1 Map of northwestern Argentina, showing the provinces of Jujuy, Salta, Tucuman, Catamarca, and Santiago del Estero 


The insert shows the relative position of the provinces within Argentina. National Parks are indicated with numbers 1 to 6 as follow: 1) 


Parque Nacional Baritú (Salta province); 2) Parque Nacional Calilegua (Jujuy province); 3) Parque Nacional El Rey (Salta province); 4) 


Parque Nacional Los Cardones (Salta province); 5) Parque Nacional Copo (Santiago del Estero province); 6) Parque Nacional Campo de 


los Alisos (Tucuman province). 
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Fig. 2 Map of northwestern Argentina showing the localities where bats and small marsupials were recorded (black points), 
as well as the phytogeographic province of the Yungas (solid gray fill) 
a) Zone of the Mountain Range of Ambato; b) Zone of the Mountain Range of Ancasti, both in Catamarca province. 


intergradation of Andean and pre-Andean mountains and 
lowland desert valleys and fragments of forests that 
protrude into temperate arid and semiarid habitats. In 
Argentina, the Yungas are bordered to the east by the 
phytogeographic province of Chaco and to the west by 
the Monte and Prepuna provinces and in some sectors by 
the Puna province (Cabrera, 1976). Altitude varies from 
300 to 3,000 m and three altitudinal districts can be 
differentiated (Meyer, 1963; Cabrera, 1976; Brown, 
1995a, b; Prado, 1995). Additionally three latitudinal 
sectors were proposed that were associated with rainfall 
and orography: North, Centre, and South (Brown & 
Ramadori, 1989; Morales et al, 1995; Morales, 1996; 
Brown et al, 2001). 
1.2 Data 

Advances in the knowledge of small mammals of 
NWA are well documented for the last several decades 


the of 
biogeographic analysis. Taxonomic arrangement follows 


allowing achievement comprehensive 
Barquez et al (2006). Taxonomic identifications were 
made by comparisons of morphology and external and 
cranial morphometrics with reference material, type 
descriptions of species, and published taxonomic 
revisions of particular species or species groups. Forty 
two species of bats have been cited for the region, 39 of 
which are in the Yungas (Barquez et al, 1993; Barquez et 
al, 1999a, b; Barquez & Diaz, 2001; Barquez, 2006; Diaz 
& Barquez, 2007). Marsupials are represented by 11 
species, nine in the Yungas (Flores et al, 2000; Diaz et al, 
2002; Flores, 2003, 2006; Jayat & Miotti, 2006; Flores et 
al, 2007). The list of species and their distributional data 
were obtained from field surveys, and from literature and 
systematic collections. One thousand, one hundred and 


eighty nine (1189) records of 52 species of NWA were 
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obtained, including the 42 species of bats and ten of the 
11 small marsupials from the five NWA provinces of 
Jujuy, Salta, Tucuman, Catamarca, and Santiago del 
Estero (Fig. 2). Five species were registered at only one 
site (Appendix 1). Some georeferences were taken in the 
field with Global Positioning System (GPS) at the time 
of collection while others were obtained from gazetteers, 
maps, or satellite images. 
1.3 Methodology 

The methodology proposed by Szumik et al (2002), 
later modified by Szumik and Goloboff (2004), was used 
to identifying the areas of endemism determined by the 
distribution of bats and small marsupials in NWA. The 
method consisted in dividing the study region into cells, 
and counting the number of species that can be 
considered endemic to a particular set of cells. The sets 
with the highest number of endemic species were 
preferred (Szumik et al, 2002; Ruggiero & Ezcurra, 2003; 
Szumik & Goloboff, 2004). The 
implemented in two computer programs, NDM/VNDM 
v.2.5 (Goloboff, 2005). 

In this study an analysis of endemicity was carried 


method was 


out by using a heuristic search. It was established that the 
sets had to be found by adding/eliminating one cell at a 
time, and that these sets had to be defined by two or 
more endemic species. Also, their endemicity value 
ought to be higher than or equal to that of two species 
with a perfect homopatrid distribution. Cell sets partially 
overlapping with more than 50 percent of species in 
common were ruled out, obviously retaining those with 
the highest endemicity value. 

Given that there are no criteria to select an optimal 
grid and that there are patterns of quite different sizes, 
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the data were analyzed using various different grids and 
database of the 
georeferenced records by using grids with cells of 


fill options. We analyzed the 
different shapes (square, horizontally, and vertically 
rectangular) and sizes (Tab. 1). The method considers 
three categories of records: absent, present (including 
inferred), and assumed (as present). A species may not 
have been recorded in a cell but its presence may be 
suspected or assumed. On the other hand, a species 
absent from one cell but present in surrounding cells can 
be inferred as present in that cell. The assumed and 
inferred records can be filled by NDM with a simple 
algorithm (Szumik & Goloboff, 2004). The program 
offers to the researchers the possibility to set the filler 
radius size as a percentage of the cell size. The radius 
sizes used for assumed and inferred presences when 
database with the different 
considered are given in Tab. 1. 


analyzing the grids 

The origin of the grid was arbitrary defined and 
located at x=31.700 and y=20.530. For the mapping of 
the resulting areas, political and phytogeographic maps 
of NWA were used. The areas of endemism and the 
record points of the species analyzed were graphed on 
the maps using the Global Mapper v8.03 program. 


2 Results 


Marsupials have been useful at characterizing areas 
of endemism, and eight out of ten species had positive 
endemism scores at different spatial scales (Tab. 3). 
Regarding bats, our results show that 23 out of the 42 
species were good indicators and permit the definition of 
areas clearly coincident with the Yungas of Argentina 
with striking accuracy (Tab. 3). 


Tab. 1 Search parameters for the analysis with grid of cells of different shapes and different sizes 


(expressed in degrees) 

















Cells Filling values 
Shape Size 00/00 40/60 60/80 80/100 100/120 120/140 100/200 

1.00x 1.00 X 
0.50x 0.50 X X X 

3 0.30x 0.30 X X X X X X X 
0.25x 0.25 X X X X X X X 
0.50x 1.00 X 

HR 0.30X 0.50 X X X 
0.25x 0.30 X X X X X X X 
1.00x 0.50 X 

VR 0.50X 0.30 X X X 
0.30x 0.25 X X X X X X X 





S: square; HR: horizontally rectangular; R V: vertically rectangular. 


The size of the radio to establish assumed and inferred presences is presented under “Filling values" (assumed / 


inferred); “00/00” indicates neither inferred nor assumed presences. The “X” indicates the combination of search 


parameters considered, this is, the different filling values used for each different grid of cells. 
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Tab. 2 Summarized results obtained in the analysis of grids with cells of different shapes and sizes (expressed in degrees) 
and considering different values of filler for assumed and inferred presences 

















Cells 
S HR VR 
shape 
Cells 1x1 0.5x 0.5 0.3x 0.3 1x 0.5 0.5x 0.3 0.5x1 0.3X 0.5 
size Score Species Score Species Score Species Score Species Score Species Score Species Score Species 
X ià 12b 3.53 4b 4.77 6b 2:32 13 b 7.08 9b 2.20 2-3 b 
l 3m 4.43 12m 1m 3.57 2m 2m 3.55 12m 
6b 2.12 12b 
B 5.69 = E = = = = = = E 
3m 2.89 2-3m 
2 is 6b 2.01 1-6b 2.05 3b 4.57 6b 2.01 2-4 b 5.74 6b 2.15 2-4 b 
í 2m 4.56 12m 2.19 1m 3.11 Im 2m 2.79 12m 





S: squared; H R: horizontally rectangular; V R: vertically rectangular. 


Under “score” we consign one value when we obtained a single area under the search parameters or the minimum and maximum values when we 


obtained more than one area equivalent under different search parameters. Under “species” we detail the number of bats species (b) and of 


marsupials species (m), and we present the minimum and maximum values when there were more than one area equivalent under different search 


parameters. A, B, C: Generalized areas of endemism. 


Several areas of endemism, remarkably coincident 
with the Yungas in almost all their extension, or with 
combinations of their latitudinal sectors, were obtained 
(Tab. 2; Fig. 3). It is worth mentioning that these results 
were almost identical, and stable, with all the search 
parameters used. The following areas were obtained: 

A) This generalized area has been treated as a 
“Peninsula” (Ojeda et al, 2008). It represents all of the 
areas of endemism obtained from the Yungas of NWA 
along their entire latitudinal extension (Fig. 3A), and it 
was obtained several times, with grids of square, 
horizontally, and vertically rectangular cells, with scores 
that range widely between 2.20 and 9.13 (Tab. 2). 

All coalescent areas of endemism are included 
within this generalized area, even if some have 
individual limits which are slightly different from each 
other. The obtained areas that can be grouped here cover 
almost all the surface of the Yungas in NWA. They 
range from 22° S latitude in the north of Salta province, 
through the east of the province of Jujuy, centre of Salta, 
and centre and west of the province of Tucuman, south to 
approximately 28.5° south latitude in the province of 
Catamarca (Fig. 3A). None of these individual coalescent 
areas includes the southern portions of the Yungas 
formation (i.e., the Mountain Ranges of Ambato and 
Ancasti, in Catamarca) (Fig. 2). 

While some of the obtained areas of endemism are 
contiguous, other have disjunctions in central Salta, 
between 24.5° and 25.5° south latitude (Fig. 3A), that 
may be related to at least two factors of importance. First, 
an important discontinuity is located between the north- 
central and south-central mountain ranges of Salta 


province, resulting in an important discontinuity of the 
forests formation in the region (Fig. 2). Second, but 
related to the first, the amount of records for the area is 
quite low if compared with the known records for other 
sectors of the Yungas in Argentina, except the 
southernmost sectors of the Yungas in the province of 
Catamarca (Fig. 2). 

The score of endemicity was given by 13 species of 
bats (Artibeus 


Eptesicus 


planirostris, Chrotopterus auritus, 


diminutus, Histiotus laephotis, Lasiurus 
blossevillii, Myotis albescens, M. keaysi, M. riparius, 
Nyctinomops laticaudatus, N. macrotis, Sturnira 
erythromos, S. lilium, and S. oporaphilum) and four 
marsupials (Micoureus constantiae, Thylamys cinderella, 
T. pallidior, and T. sponsorius) (Tab. 3, in part). 

The highest individual endemicity values for bats 
(greater than 0.67) were found in A. planirostris, C. 
auritus, L. blossevillii, M. keaysi, S. erythromos, S. lilium, 
and S. oporaphilum. Artibeus planirostris and C. auritus 
are mostly found in the Yungas in NWA, although 
Chacoan records also exist. However, these species have 
high endemicity values because the localities of the Dry 
Chaco are very close, and closely related to the Yungas. 
L. blossevillii is a widely distributed species in Argentina, 
and most of its records in the northwest correspond to the 
Yungas, with just a few in the Dry Chaco. In Argentina, 
M. keaysi and S. erythromos have been recorded only in 
the provinces of Jujuy, Salta, and Tucuman, with most of 
the records belonging to the Yungas, and just some to the 
Dry Chaco. Sturnira lilium is a very abundant species all 
over the entire north of the country and in the northwest 


inhabits primarily in the Yungas, and secondarily in the 
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Tab. 3 Characterizers species of the three generalized areas of endemism obtained in the analysis 
with grids of squared cells of different sizes (expressed in degrees) and considering different 
filling values for assumed and inferred presences (not shown) 


María L. Sandoval et al: Small Mammals Defining Areas of Endemism 








Areas A B C 

Cell size 1 0.5 0.3 1 0.5 0.3 1 0.5 0.3 
Species 

A. caudifer = a ge se = = 0.90 0.74 0.64 
A. planirostris 0.74 0.79 = 0.75 0.32 P = er = 
C. auritus 0.67 0.74 = Se = = = = = 
C. planirostris = = x 0.78 e = 0.72 0.61 = 
D. youngi = S = i = S 0.70 0.76 = 
E. diminutus 0.62 B B E = = = e = 
E. patagonicus = = E 0.18 = = = E -= 
G. soricina be zz i -= = zi 0.80 0.73 0.76 
H. laephotis 0.63 0.61 x: -= -= - = = 7 
L. blossevillii 0.44 0.70 B B s = £ B = 
M. microtis -= = =i zs = zs 0.72 = 
M. rufus = = = 0.67 = = 0.56 0.71 = 
M. keaysi 0.80 = = = x = xx = =< 
M. riparius 0.62 st 2 = = i = = 2 
N. leporinus T E ux 0.68 - P ^ i: = 
N. laticaudatus 0.47 i= 5 5 = = = E = 
N. macrotis 0.51 2 = = x = = a E 
P. bilabiatum S = = = = = 0.80 0.85 0.79 
S. erythromos 0.56 0.69 D = + = + = = 
S. lilium 0.67 0.76 = = js = = = E 
S. oporaphilum 0.67 0.70 2 0.68 0.75 = X > = 
T. bidens ES E -= = -= - -= 0.81 -= 
C. chacoensis = FE = => = = 0.70 0.79 = 
C. ignitus E = S = = = = 0.76 = 
M. constantiae = = ge 0.63 0.68 = = = = 
T. cinderella 0.84 0.85 = 0.71 0.75 = x = = 
T. pallidior 0.27 a e se = = = = a 
T. sponsorius 0.63 0.78 x E = E = = = 
T. venustus = = S = = = 0.70 0.83 = 
Thylamys sp 5 in = 0.61 = = = a = 


The endemicity values given by each species to the areas of endemism obtained are shown. In the case of a 


species that contributes to more than one equivalent area, we indicate the highest value. A, B, C: 


Generalized areas of endemism. 


67° 66° 





67° 
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Fig. 3 Map of northwestern Argentina showing the first (A), the second (B) and the third(C) generalized areas of endemism obtained 
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Dry Chaco and in the Monte Desert; other species of the 
genus, S. oporaphilum, is exclusive of the Yungas in 
Argentina. Among marsupials, M. constantiae, T. 
cinderella, and T. sponsorius have high endemicity valu- 
es, which is consistent with previous distributional 
analyses (Flores, 2006; Flores et al, 2007) that conclude 
that these species are exclusive of the Yungas. 

B) This area represents all the obtained areas of 
endemism that includes the Yungas of NWA but do not 
include their southernmost portion in the provinces of 
Tucuman and Catamarca (Fig. 3B). This area was 
obtained a few times, with grids of square cells, with 
scores that range between 2.12 and 5.69 (Tab. 2). 

As in A (the previous case) in this generalized area, 
all of the obtained areas of endemism, that occupy the 
same zone, are included. They run from the 22? south 
latitude in the north of Salta province, through the east of 
Jujuy province, and centre of Salta, to approximately 27? 
south latitude in the province of Tucuman (Fig. 3B). 
Some of obtained areas are contiguous, while others 
show disjunctions in the centre of the province of Salta, 
between 24.5? and 25.5? south latitude. 

The species that contribute to the score of 
endemicity are six bat species (A. planirostris, Cynomops 
planirostris, Eumops patagonicus, Molossus rufus, 
Noctilio leporinus, and S. oporaphilum) and three 
marsupials (M. constantiae, T. cinderella, and Thylamys 
sp.) (Tab. 3). 

In Argentina, C. planirostris has been recorded only 
in the provinces of Jujuy and Salta; its endemicity value 
is high (=0.78) because all of the confirmed records are 
from the Yungas, although this species can probably 
inhabit also the Dry Chaco (Barquez, 2006). E. 
patagonicus is widely distributed in Argentina, and have 
a very low individual value (=0.18) because in NWA its 
records are both in the Yungas as in the Dry Chaco, but 
without any predominance of records in the Yungas. 
Molossus rufus is distributed in central and northern 
Argentina; in the NWA it has been recorded mostly in 
the Yungas, with some records in the Dry Chaco. The 
records of N. leporinus in NWA are also mostly found in 
localities of the Yungas forests and very limited 
individuals are from the Dry Chaco. The only three 
known records of Thylamys sp. also belongs to the 
Yungas region. 

C) This third generalized area represents all the 
obtained areas of endemism that includes the 
northernmost portions of the Yungas in NWA (Fig. 3C), 
which are the north of Salta and the east of Jujuy 
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provinces. The retrieved areas that can be grouped in this 
category runs from the 22° south latitude in the north of 
Salta, south to the 24.5° south latitude in the province of 
Jujuy (Fig. 3C). 

The area was obtained with all the combinations of 
search parameters, with grids of square, horizontally, and 
vertically rectangular cells, and with scores that range 
between 2.01 and 5.88 (Tab. 2). The species that 
contribute to the score of endemicity are eight bats 
(Anoura caudifer, C. planirostris, Diaemus youngi, 
Glossophaga soricina, Micronycteris microtis, M. rufus, 
Pygoderma bilabiatum, and Tonatia bidens) and three 
marsupials (Cryptonanus chacoensis, C. ignitus, and 
Thylamys venustus) (partially shown in Tab. 3). Anoura 
caudifer is a species exclusive of the Yungas in 
Argentina, with records in the provinces of Jujuy and 
Salta. Diaemus youngi has been recorded in the forests of 
Misiones and in the Yungas of Jujuy and Salta in NWA. 
Glossophaga soricina has a slightly wider distribution, 
covering provinces of the northeast and centre of the 
country. In NWA all of the records correspond to the 
Yungas. Until the implementation of this analysis, 
Micronycteris microtis has been recorded in only one 
locality of the Yungas in the province of Salta. 
Pygoderma bilabiatum has been recorded in the 
provinces of Corrientes and Misiones in the northeast, 
and in Jujuy and Salta in the NWA where all the records 
correspond to the Yungas. Tonatia bidens has been 
recorded in the forests of Misiones and in the Yungas of 
Jujuy in NWA. Among marsupials, C. ignitus is known 
only by one specimen from the type locality, in the 
province of Jujuy, located in the Yungas. Thylamys 
venustus is known only from Yungas localities of the 
province of Salta. 


3 Discussion 


This study consists in a formal analysis of 
endemism based on distributional records of mammals, 
and focused on species of bats and small marsupials from 
northwestern Argentina. We decided to apply an 
optimality criterion (Szumik et al, 2002; Szumik & 
Goloboff, 2004). On the other hand, the parsimony 
analysis of endemicity (PAE) (Rosen, 1988; Morrone, 
1994) has been widely used to identify areas of 
endemism for different taxa in several regions (Escalante 
et al, 2009; and references therein). PAE uses a 
parsimony algorithm over a matrix of grid cells 
(equivalent to taxa in a classic phylogenetic analysis) by 
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taxa (equivalent to characters in a classic phylogenetic 
analysis) in order to obtain an area cladogram. In such 
cladogram, clades are considered areas of endemism. 
Despite its wide use, PAE provide a worst resolution 
than NDM and allow the identification of fewer areas of 
endemism and fewer diagnostic species (Casagranda & 
Taher, 2007; Carine et al, 2009; Escalante et al, 2009). 

Eighty percent of the species of small marsupials 
included in the analysis scored positively to the 
identified areas of endemism (eight -C. chacoensis, C. 
ignitus, Micoureus constantiae, Thylamys cinderella, T. 
pallidior, T. sponsorius, T. sp, and T. venustus- of ten 
species); this is in concordance with the fact that each 
species has a particular and restricted ecological 
distribution. Regarding bats, a priori it is usually 
considered that they may not be considered important 
indicators of endemicity, because of their ability to 
disperse by flying. Contrary to predictions, our results 
show that almost 55 percent of the species of bats gave 
score to these areas (23 out of the 42 species). Therefore, 
our results show that small mammals are efficient 
indicators of areas of endemism; they also show that 
those areas are congruent with the botanical definition of 
the Yungas of NWA and their subdivisions, an 
assumption that was not previously tested with any 
formal quantitative method; we also found that non- 
volant species are better indicators than volant at 
regional scale, but volant are better indicators than it was 
previously thought. 

According to Koopman (1982), in Argentina there 
is a massive, relatively homogeneous zoogeographic 
subregion, the Patagonian subregion, which includes 
almost all of the bats of Argentina, with an intrusion of 
the Eastern Brazilian Highlands subregion into 
northeastern Argentina. Barquez et al (1999b) have 
argued that this assumption is no longer acceptable 
because, according to their data, 1t was observed that the 
bat fauna of Argentina is more complex than suggested, 
but not sufficiently rich to delineate faunal regions. Other 
authors also believed that bats are not of great utility in 
delineating faunal regions (see Willig & Mares, 1989). In 
contrast to these studies, our results are markedly 
different and somewhat similar to what it was found in 
other regions, for example those reported by Andersen 
(Andersen, 1912) who as early as in the 1910 was able to 
determine solid biogeographic patterns on the basis of 
detailed distributions of bats in the Old World Tropics. 
Then we suggest that bats could be considered as 
important faunal elements to characterize areas of 


endemism, together with a diversity of terrestrial small 
mammals. 

In this study the areas of endemism were based on 
the points of occurrence of the species, and only records 
associated to museum voucher specimens, or highly 
included. It should be 
considered that, with respect to literature records, the 


reliable references, were 
identity of the species must be tentative. Also geographic 
biases may exist in relation to distributional points as 
they were not randomly recorded but distributed 
unevenly along roads, or in easily accessible areas. Other 
problems are the ecological biases associated with some 
aspects of the biology of the species that make it unlikely 
to be documented. In this sense, Ojeda et al (2003) and 
Tabeni et al (2004) have noted that the correct 
identification of priority areas for conservation may be 
potentially affected by the confused taxonomy of some 
groups, as well as by the unequal sampling effort in 
biomes, the obsolescence of the databases, or by changes 
of the landscape caused by human activities. This is also 
true with respect to the identification of areas of 
endemism. Fortunately, there are some alternatives to 
avoid these problems. The most commonly used 
alternatives are two: To apply the optimality criterion to 
data consisting in potential distributional models and to 
use the fill option in software NDM (see program 
documentation). Both techniques led to finding areas of 
endemism previously unidentified, and provided new 
species to diagnose the areas (Escalante et al, 2009). But 
both may fail because they may overpredict the 
distributional areas of the species. Moreover, modelling 
species distributions can involve the same biases as the 
use of points of occurrence because it makes 
assumptions about the primary data (uniform and random 
sampling, equal probability of capture of all species 
present at a site) that are not always true. Also, as noted 
by Tabeni et al (2004), using georeferenced records of 
occurrence optimizes the quality of analysis as opposed 
to the use of species distribution maps, which tend to 
overestimate species diversity. Biogeographic analyses 
based in occurrence points can be quite good despite 
existing some geographic and ecological biases in 
primary data, if a sufficiently large number of points is 
available. As Ojeda et al (2003) sustain, although maps 
of distributional areas of species are useful at larger 
spatial scales, at smaller spatial scales the analysis needs 
more detailed data. Moreover, the mammal fauna of 
NWA seems to be fairly well known for the total area 
(Ojeda et al 2003), thus allowing the implementation of 
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reliable analysis. 

In the present analysis, we have identified areas of 
endemism along the entire latitudinal extension of the 
Yungas in Argentina, as well as areas that correspond 
just to portions of the whole formation. This is in 
latitudinal 
impoverishment proposed by Ojeda & Mares (1989), and 


accordance with the pattern of 
in close correspondence with vegetation patterns. In 
particular, the Northern sector of the Yungas of NWA 
(sensu Brown & Ramadori, 1989) (Fig. 3C) was obtained 
with all the parameters (grid of cells with different 
shapes and sizes and with different filling values of 
assumed and inferred presences). Our analysis allows 
recognizing the Northern sector as the best defined sector 
because of the presence of unique species not found 
further south. This is not surprising because, as indicated 
by Ojeda & Mares (1989), the Northern sector has the 
higher species richness. Morales (1996) reported that the 
tree floras of the Central and Southern sectors are 
impoverished versions of the Northern sector. Our data 
support such impoverished pattern based on small 
marsupials and bats distribution, which is reflected in the 
nestedness of Northern areas with respect to the Yungas 
as a whole, as well as in the failure to recover the 
Southern and Central sectors of the Yungas as 
independent areas of endemism. The nested pattern of 
endemism is possible and frequent to observe. It consists 
in one or more areas contained within a larger one. This 
is the case of the Northern Yungas area of endemism. 
The species that give score to the Northern Yungas are 
different to those that give score to the larger areas of 
endemism that include it. All of them have distribution 


areas that overlap in part but do not define the same areas. 


For example, the nectar-feeding bats give score to the 
Northern Yungas area of endemism but not to the others. 
These tropical species occur here and no further to the 
south likely because of resource limitations. The 
southern limits of Glossophaga and Anoura may reflect 
the change of degree of seasonality between the forested 
fragments (Ojeda & Mares, 1989). In a similar way, 
several other species that are found only in the North 
define a small area that is the one with higher diversity 
(Brown & Ramadori, 1989; Ojeda & Mares, 1989). In an 
attempt to explain this high diversity, Brown et al (2001) 
suggested that the Bermejo River Upper Basin in 
Argentina and adjacent areas in Bolivia possibly acted as 
a biodiversity refuge during the Pleistocene. 

The main hypothesis of vicariance biogeography 
postulates that the emergence of barriers separate 
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simultaneously the distribution of several taxa, producing 
distributional patterns that can be explained as a result of 
a common history. The appearance of a distributional 
range, with regard to continuity or discontinuity, may be 
only an effect of the scale of the study. A real separation 
requires of a real barrier to the dispersal of individuals 
belonging to populations separated by it (Espinosa 
Organista et al, 2001). For example, the location of 
transitional lines among biogeographic regions, has often 
been a reason of controversy, because of the lack of well 
defined climatic, environmental, or biotic barriers 
(current or past) among regions, and/or because different 
groups of organisms used to identify the barriers, 
respond (or responded in the past) differently to the 
presence of those climatic, environmental, and biotic 
factors. In the case of the Argentine Yungas, they are 
completely surrounded by semi-arid and arid 
environments, with very different and contrasting 
climatic conditions, which could be acting as very 
important and efficient barriers. It is well known that 
abrupt changes in the fauna are coincident with similar 
topographic, climatic, and floristic changes. In the 
Yungas, the most important geographic barriers could be, 
first of all, the different environmental conditions that 
limit, in elevation, the distribution of organisms; in 
addition, and limiting their latitudinal distribution, there 
are series of geographic barriers, as the great valleys, 
between the principal mountain chains. The intermontane 
valleys, with different and often contrasting climatic 
conditions, act as barriers for the dispersion of some 
species. The transitional forests, a natural connection 
among latitudinal sectors, have practically disappeared, 
contributing to accentuate the isolation between them. In 
addition, barriers can be of various types, not only 
topographical; for example Ojeda & Mares (1989) 
proposed that low temperatures may be acting as barriers 
to the dispersal of some tropical species in the Yungas of 
northwestern Argentina. 

The establishment of protected areas requires the 
implementation of efficient methods that allow the 
selection of biodiversity hotspots, or priority areas for 
conservation in situ (Williams et al, 1996; Tabeni et al, 
2004). Hotspots are defined as those areas that 
individually contain the highest species richness, or the 
highest concentration of endemic, rare or endangered 
species (Williams et al, 1996; Tabeni et al, 2004). Mares 
(1992) pointed out that taxonomic diversity and number 
of endemic taxa are high in the Western Montane Forests, 
in which the Yungas are included. Tabeni et al (2004) 
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found that both the methods to identify hot spots through 
the analysis of richness, as well as those employing the 
analysis of endangered species, have led to the selection 
of priority areas for conservation in the Yungas. Also, 
Ojeda et al (2003) have found that the species richest 
grid cells in relation to total richness, to mega-species 
and to endemic species, were located in the Yungas. In 
fact, the most important protected areas in northwestern 
Argentina are located in the Yungas (Fig. 1). Our results 
also highlight the value of the Yungas in terms of its 
endemic species and, if the location of the protected 
areas proposed by other authors (Ojeda et al, 2003; 
Tabeni et al, 2004) is compared with ours, the 
coincidence is clear. Nevertheless, additional studies are 
needed for the Monte and Chaco biomes (Ojeda et al, 
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